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THE ORIGIN OF ^opE AND OTHER SHORT-LIVED 
RADIONUCLIDES IN THE EARLY SOLAR SYSTEM 

Matthieu Gounelle & Anders Meibom^ 



Abstract. Establishing the origin of short-lived radionuclides (SLRs) 
with half-lives < 100 Myr has important implications for the astro- 
physical context of our Sun's birth place. We review here the different 
origins proposed for the variety of SLRs present in the solar accretion 
disk 4.57 Ga ago. Special emphasis is given to an enhanced Galactic 
background origin for ''^Fe which was inherited from several supernovae 
belonging to previous episodes of star formation, rather than from a 
nearby, contemporaneous supernova. 

1 Introduction 

Short-lived radionuclei (SLRs) are radioactive isotopes with half-lives < 100 Myr. 
Their presence in the solar protoplanetary disk is inferred from excesses in their 
daughter isotopes in various meteorite components such as Calcium-, Aluminium- 
rich Inclusions (CAIs), chondrules and planetary differentiates. 

The abundance of some SLRs in the solar accretion disk is compatible with 
the expectations of continuous galactic nucleosynthesis, while some require a last 
minute origin, such as local production via irradiation in the solar accretion disk, 
or external stellar nucleosynthesis followed by injection. Additionally some SLRs 
' such as ^°Fe could be inherited from previous episodes of star formation in the 
immediate neighborhood of our nascent solar system. The origin of SLRs has many 
implications for the astrophysical context of our Sun's birth, early solar system 
chronology, stellar nucleosynthesis models or irradiation processes around young 
stellar objects. 

There have been many recent reviews covering different aspects of SLRs in the 
last years (e.g. Wadhwa et al. 2007; Goswami et al. 2005; McKeegan & Davis 
2004). In the present work, aimed at students who followed Les Houches' course, 
we wish to present the experimental situation focusing on the latest analytical 
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Fig. 1. Isochron diagram for tlie CAI MRS6 in Leoville (unpublislied data obtained 
by tlie senior autiior at UCLA in Edward Young's laboratory, using LA-ICPMS). Most 
points fall on an isochron corresponding to an initial ^^Al/^'^Al ratio of 4.6 x 10~^ 
(canonical ratio). Two spinel data points fall on an isochron corresponding to an initial 
■^''Al/^^Al ratio of 6 x 10^^ (supercanonical ratio). It is obvious from the figure that a 
large fractionation between the radionuclide and its daughter element (Al and Mg in the 
case of ^''Al) is required to establish an isochron. 



results as well as the latest advances made on the modeling front. Special emphasis 
is given to ^"Fe which plays an important role among SLRs, because it cannot 
be made by energetic particles irradiation. As such, it can help constrain the 
environment in which our solar system formed (e.g. Gounelle & Meibom 2008). 
We refer the readers to previous reviews such as Wadhwa et al. (2007) for a more 
detailed synthesis. 

2 The experimental perspective 

Short-lived radionuclei have now entirely decayed. Their past presence in the solar 
system is inferred from the excess of their daughter isotopes. For example ^^Al 
decays to ^^Mg with a half-life of 0.74 Myr. Excesses of ^^Mg correlating with the 
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aluminium abundance define an isochron and establish the past presence of ^^Al 
in the solar system (Fig. 1). 

To estimate the initial value of SLRs in the accretion disk, it is important to 
identify the most ancient phases which formed in the accretion disk. CAIs have the 
oldest Pb-Pb measured age of all extraterrestrial samples (e.g. Amelin et al. 2006). 
They are therefore believed to represent the first solids in the solar accretion disk, 
and their content in SLRs is taken as the solar system initial value. Some SLRs 
cannot be detected in CAIs, because the fractionation between the radionuclide 
and its daughter element is not strong enough to yield an isotopic effect that can 
be analytically resolved (sec Fig. 1). The abundance of these SLRs is therefore 
measured in other phases/meteorites and the abundance in CAIs is calculated 
using an other isotopic system (another SLR or the Pb-Pb age) for which data 
are available both in CAIs and in the other phases/meteorites under investigation. 
An alternative to CAIs for determining the initial solar system ratio is to use 
bulk carbonaceous chondrites (see Table 1 and text below), assuming that their 
precursors separated from the accretion disk at time zero (i.e. the CAI formation 
time) . 

It is important to mention that the initial solar system value is for some SLRs 
poorly defined (Gounelle 2006). This is because, when two isotopic systems arc 
used to calculate the abundance in CAIs (see above), it is implicitly assumed 
that the two isotopic systems record the same event. This assumption might 
be entirely flawed as, for example, diffusion coefficients can vary from element 
to element. In addition, if some SLRs were heterogeneously distributed in the 
protoplanetary disk, this very concept of an initial value would be meaningless 
(Gounelle & Russell 2005). Finally, it is worth noting that CAIs might not be 
the oldest solids, and that they might have formed contemporaneously with some 
chondrules (Markovski et al. 2006). 

In the following paragraphs, we discuss the experimental situation in detail for 
some SLRs, either because of their special interest, or because some important 
progress was made recently. For the SLRs not discussed here, we adopted the 
numbers mentioned in Wadhwa et al. (2007), and the reader is referred to that 
review for an in-depth discussion. Initial values are summarized in Table 1. Despite 
the fact that one initial value is given for each radionuclide, one should keep in 
mind that there might be a variability and that the initial value is not necessarily 
precisely known (see above). 

Evidence for the in situ decay of ^Be has been found in one CAI by Chaussidon 
et al. (2006), with an initial ratio ^Bc/''Bc « 6 x 10~^. -'^"Be has been identified in 
CAIs from CV3 and CM2 chondrite groups. Its initial abundance relative to ^Be 
varies between « 5 x 10~* and 1 x 10~^ (Liu et al. 2008 and references therein). 

The initial abundance of Al is at the center of a vivid debate;. For many years, 
it was considered that ^^Al was present in the CAI forming region at a canonical 
level (^'^Al/^^Al ~ 4.5 x 10~^, MacPherson et al. 1995). Based on high precision, 
in situ data, Yomig et al. (2005) as well as Cosarinsky et al. (2006) proposed that it 
could be closer to 6 or 7 x 10~^ (supercanonical ratio). If CAIs formed originally 
with such a supercanonical ratio, the canonical ratio would be due to resetting 
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Table 1. The initial value of SLRs in the early solar system 



R 


D 


S 


Ti/2 (Myr) 


R/S 


Notes 


^Be 


■^Li 


^Be 


<JO VJ-CX V o 


6x10""^ 


CAI 


lOBe 


lOg 


9Be 


1.5 


0.5-1x10-3 


CAI 


26^1 


26Mg 


27A1 


0.74 


5-6 xlO"'"^ 


CAI 


36C1 


36s/36Ar 


35Q 


0.3 


> 1.6x10""' 


CAI-ALT 


41 Ca 


4lK 


40Ca 


0.1 


0.1-4x10"^ 


See text 


53Mii 


53cr 


55 Mn 


3.5 


8 xlO-6 


CCs 


eopg 


60Ni 




1.5 


< exio"'' 


Irons 


92Nb 


92 Zr 


90Nb 


36 


10-5-10-3 


See text 


107pd 


107Ag 


lOSpd 


6.5 


6x10-5 


CCs 


129j 


i29Xe 


127j 


15.7 


10-4 


Chondrites 


182Hf 




180Hf 


8.9 


10-4 


CAIs 


205pb 




204pb 


15 


1-2x10-4 


Irons 


244p^ 


fission Xe 


238jj 


82 


7x10-3 


See text 



R is the radionuclide under consideration, D its daughter isotope, and S the reference 
stable isotope. CAI indicates that the initial value has been measured in a CAI. Irons 
indicates that the concerned SLR has been measured in iron meteorites and that the 
solar system initial value was calculated coupling the measured value with some other 
SLR whose initial abundance is known both in CAIs and in iron meteorites. In the case 
of ^^Cl, its abundance was measured in an alteration phase (ALT) within a CAI, and 
the solar system initial ratio was calculated assuming this phase formed 1.5 Myr after 
time zero (Lin et al. 2005) . CCs refer to initial abundances which were inferred from a 
set of carbonaceous chondrites (CCs). 



events (e.g. Young et al. 2005). A recent high precision study based on bulk CAIs 
(Jacobscn ct al. 2008) suggests that the CAIs formed with the canonical ratio 
and that the supercanonical ratio is an analytical artefact. Although the exact 
value of the initial ^^Al/^'^Al ratio matters for early solar system chronology, it is 
of little importance as far as ^^Al origin is concerned because the debate concerns 
variations of 20 %, far smaller than the precision of any reasonable model trying 
to account for the origin of SLRs. 

Srinivasan & Goswami (1994) measured an initial abimdance 4iCa/4'^Ca = 1.5 
X IQ-* in an Efremovka (CV3) CAI. However, the CAIs initial abundance of 4iCa, 
which decays into ^^K, might have been higher than previously thought. If Mg 
isotopes were reset by a secondary event, it is likely that K isotopes were too, 
lowering the now measured initial 4iCa/4°Ca ratio of CAIs. Assuming a similar 
temperature closure for K and Mg, and applying the exponential decay law to 
41 Ca and ^^Al, the initial 4iCa/4°Ca ratio would have been as high as 4 x lO-'' 
(see Gounelle et al. 2006). 

Moynier et al. (2007) recently measured the Cr isotopic composition of bulk 
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carbonaceous chondrites and inferred an initial value '''^Mn/''''''Mn = 8 x 10~^ 
for the SLR ^^Mn. This is lower than the initial value calculated by Lugmair & 
Shukolyukov (2001) based on measurements of angrites. 

The initial abundance of ^"^Fc in the solar system is not precisely known. Dif- 
ferent estimates determined by Multi Collector -Inductively Coupled Mass Spec- 
trometry (MC-ICPMS) or Secondary Ionization Mass Spectrometry (SIMS) vary 
between 5 x 10~* and 1 x 10~^ for the initial ^°Fe/'''*'Fe ratio (see the discussion 
in Gounelle & Meibom 2008). Measurements performed with MC-ICPMS have 
so far failed to detect an isochron indicative of the past presence of ^°Fe (e.g. 
Dauphas et al. 2008). The most precise work performed to date on primitive car- 
bonaceous chondrites by Regelous et al. (2008) constrained the initial ^°Fe/^^Fe 
ratio to be lower than 1 x 10~^, suggesting that SIMS work might be plagued 
with unresolved interferences. We will conservatively adopt an upper limit of 6 x 
10-"^ for the initial ^°Fe/^^Fe ratio (Dauphas et al. 2008). 

Nielsen et al. (2006) measured in iron meteorites excesses of ^^^Tl correlating 
with 204pb, indicating the past presence of the SLR ^"Spb (T1/2 = 15 Myr). The 
initial solar system abundance was calculated using I-Xe ages for iron meteorites. 
Baker et al. (2007) failed to identify an Tl-Pb isochron for carbonaceous chon- 
drites, although data from the latter fall close to the isochron defined for irons. 

3 Models for the origin of sliort-lived radionuclides 

3.1 A galactic background origin for some SLRs 

On-going nucleosynthesis by a diversity of stars (supernovae, novae, AGB stars...) 
in the Galaxy continuously replenishes the interstellar medium with freshly made 
SLRs. At a given time in the history of the Galaxy, the background abundance 
of a given short-lived radionuclide R relative to its stable reference isotope S will 
depend on a diversity of parameters such as the number and nature of nucleosyn- 
thetic events responsible for the production of R over the last few half-lives of 
R, the number and nature of nucleosynthetic events responsible for the produc- 
tion of S during the history of the Galaxy, the respective yields of R and S in 
the nucleosynthetic events aforementioned, astration, the mixing timescales and 
processes of the different phases of the interstcillar medium, the rate of decay of 
R... Final isolation of the average interstellar medium from nucleosynthetic events 
introduces an extra parameter, the isolation time A, during which R decays with- 
out further addition of freshly made matter. Because two different SLRs might 
originate from different nucleosynthetic events, A varies from radionuclide to ra- 
dionuclide. In addition to the number of poorly constrained parameters mentioned 
above, a further complication occurs for those SLRs whose half-life is smaller than, 
or comparable to, typical recurrence times of relevant nucleosynthetic events. In 
such a case, granularity of nucleosynthesis becomes a key phenomenon (Meyer 
& Clayton 2000), and it becomes difficult to model a steady-state background 
interstellar medium abundance. 

Despite the uncertainties discussed above, Wasserburg et al. (2006) present 



6 



Title : will be set by the publisher 




Fig. 2. Abundances of SLRs injected by a nearby SN in the nascent solar system com- 
pared to the value observed in meteorites. Because the injected abundances depend only 
on two parameters (the mixing fraction / and the decay time A, see Eq. (1) of Gounelle 
& Meibom 2008), the calculation assumes that ^''Fe and ^'^Ca are delivered at the cor- 
rect abundance. Yields of massive stars are from Rauscher et al. (2002). Different lines 
correspond to supernovae of different masses. 



some estimates of the background interstellar medium abundance of SLRs. Ac- 
cording to these authors, ^^''Pu and ^*^Hf had solar system abundances compat- 
ible with the background interstellar medium value, while ^''Be, ^^Al, "^^Cl, ''"'^Ca 
and ^"Fe were in excess in the solar system relative to the expected background 
interstellar medium abundance, calling for a last minute origin. They use the 
abundance of ^^^I, which is a r-process only nuclide, to calculate the isolation time 
(A ^ 70 Ma) since the last r-process nucleosynthetic event. If the same time 
delay is applied to ^"^Pd, which is also a r-process nuclide, they find that ^"^Pd 
is overabundant in the solar system, and requires a last minute origin, which is 
attributed to the contribution of an AGB star by Wasserburg et al. (2006). 

The estimated solar system abundance of ^°^Pd and ^^^I can however be recon- 
ciled with expectations from the continuous galactic nucleosynthesis. Still assum- 
ing that ^"^Pd and ^^^I are made in the same nucleosynthetic site, and considering 
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an isolation time of A = 43 Ma, we calculated, using the Wasserburg et al. (2006) 
numbers, that ^"^Pd was underabundant in the solar system by a factor of 3, while 
^^^I was overabundant in the solar system by a factor of 3. Given the uncertain- 
ties of the models and of the initial abundances of SLRs, it is quite reasonable 
to conclude that both ^"^Pd and "'^^^I have abundances in line with that of the 
steady-state interstellar medium background. 

Nielsen et al. (2006) and Wadhwa et al. (2007) proposed that ^ospb and 
^^Nb also originated from the galactic background. It could also be the case of 
53Mn (Wadhwa et al. 2007). Berylhum-7, i°Be, ^^k\, ^^G\, ^^Ca and ^Ope cannot 
be explained by the continuous Galactic nucleosynthesis, and therefore require a 
special origin. 

3.2 An enhanced galactic origin for ^^Fe: Tiie SPACE model 

Iron-60 cannot originate by irradiation (Lee et al. 1998). As it is unlikely a nearby 
SN delivered it in the nascent solar system (see §3.3), it has recently been proposed 
that it is inherited from previous episodes of star formation (Gounelle and Meibom 
2008). 

In the last decade, it has been proposed that molecular clouds (MCs) are 
transient, dynamically evolving, dense ISM features produced by compressive 
motions of either gravitational or turbulent origin, or some combination thereof 
(Ballesteros-Paredes et al. 2007). MCs contain a mixture of atomic and molecular 
gas (e.g. Goldsmith et al. 2008) with density ranging from 1 cm^"^ to 10^ cm^'^ 
with an average value of ~ 100 cm^'^. This MC dumpiness means that in Giant 
Molecular Clouds (GMCs), some regions might actively form stars, while others 
are diffuse and undergo a latency phase (Elmegreen 2007). 

OB associations are the outcomes of the star-forming process in molecular 
clouds. OB associations are divided into subgroups of different age (Blaauw 1964). 
A famous example is the Scorpio-Centaurus region which consists of the Lower 
Centaurus Crux (LCC, ^ 16 Myr), the Upper Centaurus Lupus (UCL, ~ 17 Myr) 
and the Upper Scorpius (Upper Sco, 5 Myr) subregions (Preibisch & Zinnecker 
2007; Mamajek et al. 2002). The range in ages within a given region has long 
been interpreted as indicative of sequential star formation (e.g. Elmegreen & Lada 
1977). 

Relatively high concentrations of ^''Pe and other radioactivities with half- lives 
> 1 Myr are expected in second generation star-forming regions. This is because 
supernovae explosions from an older region, e.g. the LCC, enrich a younger region, 
e.g. the Upper Sco, with their nucleosynthetic products. It is therefore suggested 
that ^"^Pe in the solar system was inherited from previous episodes of star formation 
(SPACE modeQ, Gounelle et al. 2008). 

Two different astrophysical settings for the SPACE model are envisioned 
(Gounelle et al. 2008). In the first setting, it is considered that the entire molecu- 
lar cloud is assembled via turbulent convergent flows (e.g. Hartmann et al. 2001). 
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In this scenario, the gas which will make the bulk of the MC is swept up by winds 
from massive stars and supernovae explosions. In an alternative setting, we assume 
that supernovae explode in a pre-existing molecular cloud which is large enough 
(such as a GMC) to have regions evolving at different paces. Winds from massive 
stars belonging to an older region accumulate a dense shell of gas as in the collect 
& collapse model (Elmegreen & Lada 1977). When massive stars explode as SNe, 
star formation is promoted in the dense shell. In both cases, ^"Fe present in the 
SNe ejecta is delivered while dense gas is accumulating. Because star formation 
occurs rapidly (on a Myr timescale) once the gas is dense enough, ^°Fe has 
no time to decay and is expected to be alive in the newly formed protoplanetary 
disks. 



3.3 Nearby supernovae models 

Of models based on a supernova (SN) origin for SLRs there are two types. Either 
a SN injects freshly synthesized SLRs into a nearby molecular cloud core, trig- 
gering its gravitational collapse (e.g. Cameron & Truran 1977), or directly into a 
nearby protoplanetary disk (Chevalier 2000; Hester & Desch 2005; Ouellette ct al. 
2005). The first SN scenario is now considered less likely because only very spe- 
cific conditions allow a supernova Shockwave to trigger the collapse of a molecular 
cloud core and, at the same time, inject SLRs (e.g. Boss & Vanhala 2000). In the 
second scenario, which is currently receiving a lot of attention (Chevalier 2000; 
Hester & Desch 2005; Ouellette et al. 2005; Ouellette, Desch & Hester 2007), the 
SN has to be very close (^ 0.3 pc) to the protoplanetary disk in order to allow 
the disk to intercept enough SN ejecta to account for the solar system inventory 
of SLRs. It is thus assumed that the massive star, which evolved into a SN, and 
the protoplanetary disk were coeval and formed in the same stellar cluster (e.g. 
Hester & Desch 2005). 

The probability of direct injection of SN materials into a nearby protoplanetary 
disk has recently been estimated to be less than one in thousand (Williams & 
Gaidos 2007; Gounelle & Meibom 2008). A similar low probability can be assigned 
to the molecular cloud core model. In both cases, the evolution timescale of 
massive stars is too long compared to low-mass stars evolution timescales, given 
that SLRs incorporation happened in the earliest phases of the solar system. UV 
radiation emitted by massive stars create an HII region which can extend to a few 
pc or more after a few Myr of evolution (e.g. Reach et al. 2004), preventing star 
formation in that region. Low-mass stars formed in the enrichment zone (< 1.6 
pc) before the clearing of the molecular gas would have formed planets by the time 
of the first SN explosion (^ 5 Myr after the onset of star formation). Low-mass 
stars formed around the time of the SN explosion would be too far away from the 
massive stars (at least a few pc, more likely 10 pc) to have their disks or cores 
contaminated by ^"Fe at the level observed in the solar system. 

Regardless of its astrophysical implausibility, the injection of SLRs by a nearby 
SN also fails to quantitatively satisfy the cosmochemistry data (Fig. 3). 
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3.4 Irradiation models 

Irradiation of nebular gas and/or dust by accelerated hydrogen or helium nuclei 
can result in the production of short-lived radionuclides. ^°Be has an irradiation 
origin since it cannot be made in stars (e.g. McKeegan, Chaussidon & Robert 
2000). Beryllium-7 which has been identified in one AUende CAI (Chaussidon, 
Robert & McKeegan 2005) is an unambiguous tracer that some irradiation took 
place in the early solar system. Ubiquitous, thousand-fold enhanced, and flare-like 
X-ray activity of protostars provides firm evidence for the existence of accelerated 
particles in the vicinity of the early sun (Wolk et al. 2005). The issue is whether 
some other SLRs, such as ^^Al, are co-produced during irradiation events. Some 
models, dedicated to calculate irradiation yields of SLRs, established that it is 
possible to produce "^Be, ^°Be, ^^Al, ^^Cl, ^^Ca, ^^Mn at abundances in line with 
that of the early solar system (Leya et al. 2003; Gounelle et al. 2006) provided 
that proton and helium nuclei are accelerated by impulsive events. Assuming that 
^^Al was ubiquitous in the entire protoplanetary disk, Duprat & Tatischeff (2007) 
proposed that irradiation could not account for the canonical ^^Al/^^Al ratio. 
There is however no positive evidence for ^^Al ubiquity in the protoplanetary 
disk. 

4 Conclusions 

Some SLRs such as ^^Mn, ^^Nh, ^o^Pd, ^^^m, ^o^Pb and ^'^^Pn might origi- 
nate from continuous Galactic nucleosynthesis (Wasserburg et al. 2006; Meyer & 
Clayton 2000). These correspond to the SLRs with the longest half- life. SLRs with 
intermediate half-life (^^ Al, ^°Fe) might be inherited from previous episodes of star 
formation in the molecular cloud progenitor from our solar system (Gounelle et 
al. 2008). '^Be, i°Be, ^^Cl, ^iCa might have an irradiation origin. 
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